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INTRODUCTION 


Th«  MtheMtical  description  of  electrocheaicsl  processus  at 
stetionery  electrodes  yields  psrtlsl  differentisl  equstions  of 
the  type 


for  eseh  species*  where  f(ct  is  detereined  by  the  hoeoqeneous 
kinetics  of  the  electrochesicslly  produced  intereediatea  and  the 
substrsbe  snd  i  is  s  distance  coordinate.  The  initial  condition 
notaally  is  that  the  substrate  is  the  only  species  present  in 
solution  and  is  homogeneously  distributed  in  the  k  direction. 

Th*  boundary  conditiofu  are  dictated  by  the  experiaent 
perforasd  and  by  the  kind  of  electron  transfer*  i.e.  reversible 
or  irreversible,  it  has  been  shown  for  a  variety  of  boundary 
conditions  (i«4]  that  orthogonal  collocation  (S«6|  provides  a 
powerful*  fast,  and  easily  applied  method  to  simulate  kinetic 
diffusion  problems  of  this  kind,  up  to  now*  however*  only 
boundary  conditions  related  to  reversible  electron  transfer 
reactions  have  been  used. 

We  now  report  the  theory  for  application  of  this 
mathematical  method  to  quasireverslble  electrode  mechsnismsi 
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where  the  electrode  reaction  of  A  to  b  is  described  by  a  standard 
potential  £**  a  standard  rate  constant  for  electron  transfer  k^* 
the  transfer  coefficient  a  and  the  number  n  of  transferred 
electrons,  product  B  may  reset  in  a  chemical  reaction  13)  with  a 
homogeneous  rate  constant  k. 

INITIAL  AND  BOUNDANy  CCMOITIONS 


To  simulate  a  reaction  scheme  12)-(3)  ws  have  to  solve 
equations  of  the  type  (1)  for  three  species  A*  B  and  C  (assuming 
equal  diffusion  coefficients  for  each  species): 


D  ^ 

at  ■ 

T>  ^ 

at 

at  ■ 

^  ax* 

(  «) 

(  S) 


(  C) 


The  calculation  of  Cq  is  not  essential  to  simulate  purely 
electrochemical  respu^tses  Cor  this  reaction  scheme.  However*  if 
one  wishes  to  do  so*  this  calculation  may  be  included  without  any 
problem. 

The  initial  conditions  are: 

00  : 


B 


C  A 


(  I) 


(  7) 


th»  boundary  condition*  at  an  Inflntta  dlatanca  Iron  tha 

olactrodt  arc  ^Ivan  by 

t>«, 


uharcaa  at  the  alactrode 

i>p,  kaO: 


/)tnO 

bold*,  fiara  |7.f| 


•  *J-  evf  t(4-«)»TAT  (F-F*)} 


Ull 


azt 


un 


4^«  4^- top^-cokT/nT  (F-F»)] 


U4I 


introducing  Nldiolaon**  daflnltlons  |8»91  Cor  potential  acan 
condition*  (cyclic  coltaMatry) 


,-«L* 


s.ct)^ 


lUb) 


atT 


(uhara  %tart  atartinq  potential  of  tha  potantial  auaop 

ulth  awaap  apaad  a  and  t  ^  ia  tha  tin*  aban  tha  awaap  la 
r*var*adl»  aquation  (121  can  ba  chanqad  to 


•m 

bacauaa  apacla*  C  la  not  alaetroaetlaa»  It*  flua  at  a  *  9  la 
aluaya  aarot 

(9ee  1 

Tba  ayatan  of  partial  dlffarantlal  aqvatlona  (ai-li)  Ha*  to 
ba  aolvad  undar  tba  initial  condition*  (71  and  (91  and  tha 
boundary  condition*  (91»(lll»  (191  and  (19). 

TbMISrOMIATlON  INTO  h  DINCNSIONLBSS  tftTW  09  BQUHTIONO 

NC  Hava  to  tranafow  tbla  ayatan  of  partial  dlffarantlal 
aquation*  and  tba  condition*  under  ubicb  It  la  to  aoWa  Into  a 


i 

1 


6 


diacnaionles*  tor«.  Wt:  uite  the  followirhj  transforMtion 
aquations: 


T  .«.* 

(20) 

X'x/t 

(21) 

fA*“=AAA- 

(22) 

C»*  V‘a* 

(23) 

cJ  *c<Aa* 

(24) 

3i 

12S) 

y  *  Aj/yoiy 

(24) 

w  V/m.L* 

(27) 

wh«r«  L  U  a  distanca  (roa  tha  aiactroda  uhara  no  diffusion 
occurf  dUirinq  tha  aiaulation.  thus«  tha  diaunsioniaas  systaa  is 


JT* 


9ci 

3t‘ 

ii. 

ar 


A: 

ax* 

ax* 

?V 

ax* 


(21) 

<  2$) 

OO) 


undar  tha  conditions 

T'.o,o«X<l: 


r'xj,  x^-t' 


ci--c*>o 


T*>0,  Xeo; 


*  o 

Cc^Coj')-  c,»(o,T*)S;^tT9eJ 


with 


r  CitpC-T*)  T'<Tj[ 

l«rp(T'-jr;)  Ti<T**aT; 


U3) 

(  34) 

(  3Sa) 

(  3$b) 

(34) 


(37a) 

(37bl 


DX5C.**et)tAttGN  WITH  OHIHOGOMAl.  COLLOCATION 


NOW  ws  usa  orthoqoaaX  collocation  to  dorian  a  systaa  of 
siaultsnaouB  ordinary  diCfsrsntial  aquations  froa  aach  of  tha 
threa  partial  diffsrsntial  aquations  (2t)<>(10).  Ms  fit  tha 
coefficients  of  a  polynoalal  in  c*(Sit)  such  that  tha 
diffarantial  aquations  are  fulfilled  ssactly  at  certain  points. 

It  has  been  shown  that  a  very  officiant  collocation  asthod  can  ba 
constructed  if  wa  usa  orthooonal  polynoaials  and  if  those  points 
are  tha  saroas  of  tha  polynooial  1A|.  Purthsraore.  in  this  case* 
tha  distance  derivatives  of  tha  concentration  at  these 
*collocatian  points*  are  given  by  (i.S.il 


c5.c;*c? 


(31) 

(32) 


(31) 


dl*c*  I  ^ 

tW,.' 


where  the  A|^|  and  Si^j  are  eatrls  eleecnta  only  dependent  on  the 
kind  and  decree  of  the  puiynoalal  choeen.  h  fnlynoetal  of  deqree 
N  haa  N  roots  in  the  open  Interval  between  eero  and  one  and  we 
have  to  include  the  boundary  points  at  s  •  0  and  «  •  1.  Thus» 
the  duMny  variable  j  varies  free  l  to  N*]. 

Substituting  the  second  darivatives  with  respect  to  a  In 
equation  f2t>  we  qet 
J-u  i  •^Jt 


Jr*  i  * 

*1  i*4 


Mmc*  tiH  CM  te  MpaaM  •M  tl«  bounlarr  emidltlon*  Ul)  ^ 
Introdvndi 

Jitj  I  r  . 

'•  t*4  •* 

tlw  eoncantcatian  of  apaclaa  a  at  ttia  alactroda  auctaea  la  oot 
knom  aaplleitly  but  aoy  ba  datlaad  froa  tba  boundary  eondltlona 
Ilf)  and  ualny  I  ]f ) ,  uquatlun  lit)  ylaldai 

"al  oaj 

fsa 

“  -  -  Wed 

area  «**  <131 


equation  ( 3Si  qlves  faqaln  ualnq  |JS)| 

91*1  t 

~  S;,iW’‘^e'^Csf  (»rT*j-(;(v'is>fT'>e3 ,,,, 

or  after  eapannlon  of  the  aun*  collection  of  the  terae  eontelnlnq 
c)^tO»T*|  and  solution  for  c^fdyTMt 

equations  U3)  and  can  be  solved  for  CAtObTMt 

^  AvCAm* 

w#a  a*er 

Ca^*  z\j  v’V  <«. 

bow  it  la  possible  to  subatltete  tor  ca(*»Y*>  i**  equation  141) 


and  ue  finally  qet 


”  cr«# 

[ov  z  A^  T*n-^ 

r*  a  ItTI 

Z*sj«*  ••'iJ’J  I 

It  follows  fron  equation  ( 2S>  after  insertion  of  iiq)t 
eipansion  of  the  sun  end  Introduction  of  the  boundary  condition 


'>•  >*l  J 

Again,  cs(0,T*)  any  b*  calcuUtad.  li  m  inaart  (43)  into  (4S) 
tM  arriv*  at  ^ 

[t*vi~V''*''**’®'****3  (491 

l*> 

«.  Oi  My  b*  aubitltutad  for  c2(0<T‘)  In  aquation  (4»,  ylalOlnq 


tha  tiaa  darivativa  of  tha  ooncantr«tion  of  apaciaa  B  at  cha 
collocation  pointat 

^l.A 

r  B 


collocation 

•»M  •  ^  _ 

Cwg-'M'k^njJ^  Z Vi  cf (lijo) ‘  Ti c*(,;f) 

TO  aspraas  tha  cocraaponding  darivativa  foi  apaciaa  c  m  inaart 
(39)  into  aquation  ()0I«  aapand  and  uaa  tha  boundary  condition 
I  34)} 

J  •  f  baa 


Tha  eoncantration  of  apaciaa  c  at  tha  alactroda  aurfaok.  (a  •  01 
•ay  ba  darivad  froa  aquation  (39b)  aa 

Ss>.r).- —  c;(v  T') 

M  I’sa 


aiki  uuad  to  aapraaa  tba  tiaa  dapandanca  of  cc  finally} 

dlT  I  ^a}  >*4  • 

squatlona  (47),  ($0)  and  (S3)  daflna  a  ayataa  of  3  ■  N 
ordinary  diffarantial  aquationa  uhich  nay  ba  aolvad  by  eooputar 
aubroutinaa.  Ha  nota  tha  foroal  ainilarity  of  thaaa  aapraaaioAO 
with  tha  corraapondlng  aquationa  darivad  for  an  BC  nacbanian  with 
a  ravaraibla  alactron  tranafar  raaction  (2)  undar  cyclic 
volta«aatric  conditiona.  Tha  aquationa  darivad  hara  Cor 
quaairavaraibia  alactron  tranafara  nay  ba  uMd  alao  for  tha 
ajotfiation  of  aera  coopiicatad  alactroda  raactlona  with  aultiatap 
hoooqanaoua  raactlona  or  raactlona  of  tha  BCB  typa  by  chanqing 
tha  hinatic  tar«  li.a.  f(c)  in  aquation  (1)1  to  tha  approprlata 
aapraaaion.  Thia  oakaa  tha  application  of  orthogonal  collocation 
aapaeially  aaay  onca  tha  baaic  aquationa  for  a  apacific  boundary 
valua  problan  hava  baan  darivad  and  ia  ona  advantaga  of  thia 
■athod. 

CAICULATION  OT  THK  CURABNT 


Tha  current  through  tlta  working  elaetroda  ia  qivan  by 

LzkVA3>(^\ 


fjia^O 


(94) 


Ha  haw  ahotM  (2),  that  Nicholacn  and  Shain*a  currant  function 
x(at)(9|  can  ba  calatad  to  the  dinanaionleaa  quantitiaa  «iad  «n 


tut  trariiifornttio»i  t^wtiow  I30>*(241  by 


«T  AtJVJJJ' 


Tbtt#  ctlctlatiob  of  1  or  ifat)  rtquirtt  the  celeulatloii  of 
the  first  dsrlvatlvs  of  ej^  with  rtapeet  to  X  at  the  electrode 
sorfact.  Ibis  oay  be  dona  Croo  t^tlon  <141  ualnq  equations 
U«)  and  f4«l.  It  (ollom 

C'jbvwn^cA  Vi4l* 

150 

orrtntSATiON  or  ym  Dtiinisio(ittf*a  MRMmii  • 

la  the  ca«t  of  a  roeerslblt  eloetrodt  reaction  it  hat  been 
stKam  that  the  ealaa  of  the  dtoansloniaas  paraoattr  |  hat  a 
alqnttleant  Inflotnet  on  tht  accuracy  of  tht  alaulatton  (31. 
alto*  a  ntthod  to  optlnlse  thto  value  In  terms  of  stability  of 
the  solution  of  tht  dtffarantlal  aquation  systan  lias  bean  qtven 
1101 . 

Startinq  with  nilston't  atabilUy  criterion  |ll|  for 
wa— lwq*a  pradlctor-corractor-lntaqratlon  oethod 


o.ts 


IM 


uhlch  ralatao  tht  msIo 
absolata  value  of 


stapwldth  of  the  Inteqratlon  h  to  the 


3(f) 


(SO) 


for  equatlona  of  type  (20-10)*  it  la  possible  to  calculate  a 
nasinun  value  of  0  without  caualnq  InttablJIty  of  the  solution, 
inequality  (S7>  oust  be  fulfilled  for  all  spaeies  and  collocation 
points,  frnn  equation  ( ST)  follows 

AvCv.- 3  * 

In  analoqy  equation  (SO)  yields 


While  equation  1$))  oives 


(00) 


(Oil 


As  In  the  case  of  a  reversible  electron  transfer*  ut  have  to 
nalte  sure  that  (S7)  is  fulfilled  at  t)>e  eatreoe  points  of  a 
possible  scan,  l,e.  S^(T*)  •  0  (T*  ■  T*^>>0)  and  S|(T*I  •  I  (T* 


Inequality  (S7)  would  yield  inequalities  which  ace  difficult  to 
solve  for  s*  If  we  assune*  however*  that  the  scan  starts  at  a 
potential  tmll  before  the  standard  potential  (which  is  the 
qeneral  eiperinental  procedure)  ue  can  set 


u 


IS 


0»  4 


Thus*  for  a  >  0 


0^»‘l 


(62) 


(6ia) 

(6ib) 


(64) 


(6S) 


equ«tlon«  (59)  *r^  (60)  reduce  for  S^<TM  •  0  to 

end 

)Y'A(-^A,itB.;i}.7i  Co,..yWt-. 

whlU  tor  Sj(TM  •  1 

^!7  <«i 

end 

^>■^1"  A^.v^p*  Av*Bi,;]->!  i^2,.  <«7» 


finelly  stebility  (e  neinteLoed  if  inesuelttiea  (68)-(70) 
ere  fulfijjed 


065" 


(60) 

OfcS- 

(691 

0.6  S’ 

(70) 

pointed  out  for  the  reversible 

case  tlO) «  9  hes 

to  be 

cboeen  ee  ler^e  ee  poeeibLe  to  9ive  reliebte  rceulte* 
KBSULTS  OP  niK  CALCULATIOMS 


Celculetione  uein^  the  theory  derived  Above  heve  been 
perforaed  on  en  Madehl  oa/V6  coi^Yuter  veln^  e  dleloq  pro^ree  lor 
dete  input**  funetione  heve  been  Le9endro  polynoeieie  vith 

different  degree  N»  Tebie  1  elMMe  the  vetietion  of  peek 
potentiei  seperetlon  vith  n  coeipered  to  the  *refereAce*  velue 
given  by  hicholeon  It)*  The  integretion  step  width  wee  0*5  (eV). 
i.e.  the  eieuieted  eyelio  volteeeogrea  between  O.t  end  0.3  (V) 
wee  conetructed  fro*  1200  dete  points,  the  cpu-tiee  needed  for 
different  n  i»  elao  given  in  Table  i. 

Per  b  >  9  the  peeh  potentiei  seperetlon  is  essentielly  the 
Aeee  es  given  by  btcholson*  the  difference  being  0.5  (eV)  *  l.e. 
the  sene  velue  ee  the  integretlon  step  width.  Por  eaeiler  N  eCp 


copies  of  the  progree  which  elao  concetna  other  electrochenicei 
reaction  nodela  ere  eveilebl#  oft  regueat. 


t 

1«  \ 

I 

I 

d««Ut««  considerably.  Tbe  cpu-tiee.  however,  decreases.  \ 

hccordtn9  to  this  result,  all  further  slauletlons  have  been  | 

perfonMd  with  N  •  9.  Table  2  shows  variation  of  &Cp  wKh  the  | 

klnet*e  paraaeter  V*  which  Is  connected  to  Nicholson's  | 

paraaeter  t  IBI  by  j 

nu  1 

! 

the  aaslaua  deviation  fros  Nicholson's  data  la  1.9  (aVI.  a  i 

very  low  value  conslderinq  that  AB,  the  integration  step  width,  < 

is  O.S  (aVI .  ft  should  be  noted  here  that  Nicholson's  values  are  | 

siaulated  as  well.  Nith  spprosiaataly  2*9  cpu>seconds  needed  for  . 

1290  dats  points,  tha  proqraa  works  fast  and  efficient.  I 

CONCLUSION  * 

The  technlcfue  of  orthooonal  collocation  can  be  applied  to 
sUttlats  quastreveraible  electrode  processea  under  potential 
sweep  conditions.  The  equations  dsscrtblnq  the  tiae  dependence 
of  the  eoneentraklon  of  all  apecles  are  fonMlly  very  slailar  to  | 

those  derived  (or  reveraible  electron  transfere.  The  ' 

dlaenslonlesa  paraaeter  •  can  be  optiateed  In  teraa  of  stability  I 

of  the  aolution  of  the  differential  equation  system.  The  t 

nuiserical  results  calculstsd  for  tha  *pute*  quaslreversible  i 

aleetron  transfar  (l.e.  without  hoaoqeneout  kinetic  I 

coaplieationa)  ace  very  close  to  the  reference  values  qlven  by 
Nicholson  IBI  and  coaputed  by  a  totally  different  aethod,  i 

The  application  of  this  work  Is  obviously  the  Interpretation 
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of  aore  coapHcated  electrode  reaction  aechanlsas  where 
hoaoqeneous  chealcal  reactions  are  coupled  to  the  quaslreversible 
electron  transfer.  The  eapanslcn  of  the  theory  qlven  here  to 
those  cases  Is,  as  usual  In  orthoqonal  collocation  calculations, 
easily  done.  Sinulatlon  and  coaparlson  to  eiperlaents  with  real 
substrates  underrjolnq  coapHcated  reaction  sequences  will  be 
shown  In  a  further  paper. 


\  • 
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orthoqonal  collocation,  variation  of  f*. 
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